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ABSTRACT 

We report the results of radiation-magneto-hydrodynamics calculations in 
the context of high mass star formation, using for the first time a self-consistent 
model for photon emission (i.e. via thermal emission and in radiative shocks) 
and with the high resolution necessary to resolve properly magnetic braking ef- 
fects and radiative shocks on scales < 100 AU. We investigate the combined 
effects of magnetic field, turbulence, and radiative transfer on the early phases 
of the collapse and the fragmentation of massive dense cores. We identify a new 
mechanism that inhibits initial fragmentation of massive dense cores, where mag- 
netic field and radiative transfer interplay. We show that this interplay becomes 
stronger as the magnetic field strength increases. Magnetic braking is transport- 
ing angular momentum outwards and is lowering the rotational support and is 
thus increasing the infall velocity. This enhances the radiative feedback owing to 
the accretion shock on the first core. We speculate that highly magnetized mas- 
sive dense cores are good candidates for isolated massive star formation, while 
moderately magnetized massive dense cores are more appropriate to form OB 
associations or small star clusters. 

Subject headings: magnetohydrodynamics (MHD) - radiative transfer - methods: 
numerical - stars: formation, kinematics and dynamics, massive 
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Introduction 



Massive star formation (M* > 8 M Q ) is one of the most challenging astrophysical 
problems. It is established that most massive stars form from massive prestellar cores 
and occur in high-order multiple systems (e.g. IZinnecker fc Yorkd 120071 ). Nevertheless, 



all theoretical numerical models to date show that massive prestellar cores are unlikely 



to form w ithout first fragmenting i nto several objec ts. 



work (e.g. 



Hennemann et al. 



2009 



Bontemps et al. 



2010 



n addition, recent obse rvational 



Longmore et al. 



201ll ) suggests 



that collapsing massive dense cores are less fragmented than what numerical calculations 
produce although the limited observational resolut ion available precludes a definitive 
answ er. For instance, isother mal simulations (e.g. 



2010 



Bonnell et al 



2006) . radiative calculations 



Hennebelle et al. 



2011 



Krumholz et al. 



Peters et al. 



including radiative ionisation ( iPeters et al. 



2007. 



2011 



2001 



Bonnell fc Bate 



2009 ), magneti zed ones (IWang et al. 



Seifried et al. 



20111 ) and even calculations 



20101 ). tend all to form several fragments. 



Indeed, both radiation and magnetic fields tend to reduce the number of fragments (e.g. a 
factor of about 2 for highly magn etized cores) without suppr essing the fragmentation. To 



balance this fragmentation issue, 



Krumholz fc McKee 



( 120081 ) proposed a column density 



threshold for massive star formation. However, this model only applies to massive star 
formation under certain conditions, and in particular to massive stars, whose formation has 
been possible thanks to the radiative feedback of low-mass protostars. 

In this Letter, we perform full radiation-magneto-hydrodynamics (RMHD) calculations 
of massive, turbulent, and ma gnetized dense cor e s. Th is work is an extension towards 



Commercon et al. 



(120101 ). who investigated low- mass star 



higher masses of the study by 
formation. 

The paper is organized as follows: In section 1 we discuss our numerical method and 
initial conditions. Our results are presented in section 2. In section 3 we dicsuss a new 
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scenario for massive star formation and section 4 concludes the paper. 



2. Model 
2.1. Initial conditions 



Hennebelle et al 



(1201 lh . except that 



Our initial setup is identical to the one used in 
we do not use a barotropic equation of state. We consider 100 M spherical dense cores, 
which are threaded by a magnetic field parallel to the x-direction. The initial radius of 
the sphere is 0.67 pc and the total box-length is 2.76 pc. The density profile is given by 
p(r) = p c /(l + (r/r ) 2 ), where p c = 1.4 x 10~ 20 g cm -3 is the central density, and r ~ 0.22 
pc is the extent of the central plateau. We impose a density contrast between the central 
density and the edge density of 10. The initial temperature of the core is uniform and 
equals 10 K. Outside the cloud, the matter is also set to 10 K. The adiabati c index is set to 



" y = 7 /5. We use the Rosseland /% and Planck Kp mean opacities derived in 



Semenov et al. 



( 120031 ). At temperature > 1000 K, we impose Kp = = 0.01 cm 2 g _1 , in order to limit the 
grain evaporation effect and to account for an inertia of the evaporation. We impose an 
initial subsonic turbulent velocity dispersion which follows a Kolmogorov power spectrum 
P(k) oc k~ 5 ^ 3 , where the phases are randomly sorted in the Fourier space. Only one 
realization is explored in this study. The turbulence is not artificially sustained but, does 
not really decay as the gravitational time is typically shorter than the crossing time. The 
kinetic energy power spectrum peak roughly corresponds to the box size. The ratio of the 
turbulent to gravitational energies is given by ccturb- We do not include explicitly rotation 
but the turbulent field contains angular momentum (both local and global). Note that our 



initial conditions exhi 



( iGirichidis et al 



jits a relatively flat density profile, which should favor fragmentation 



20111 ) 
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The magnetic intensity is set by the parameter \x = ( M /<&) / (M /<&) cr \t, which 



the value of mass-to-fl ux over critical mass-to-flux ratio flMouschovias fc Spitzer 



repres ents 



1976 



Hennebelle et al. 



20111 ) . In this Letter, we explore 3 magnetization degrees: \i = 130 



corresponding to a weak initial magnetic field, and \x = 5 and /i = 2 which are close to 



the observed values (e.g. 



Falgarone et al 



20081 ). We also investigate an almost spherical 



case with no magnetic field and no turbulence (model SPHYDRO) to serve as simple 
reference with respect to which the other simulations can be compared. All the simulations 
parameters are summarized in Table [TJ Future work will imply further investigations on the 
effects of turbulence or rotation, but this goes beyond the scope of this Letter. 



2.2. Numerical method 



We use the adaptive mesh refinement code RAMSES (ITevssierl 120021). which integrate s 



the s elf-consistent equations of RMHD usi ng ideal MHD (IFromang et al 



20061 ). and flux limited diffusion for RHD (jCommercon et al. 



2011a). 



2006 



Teyssier 



We impose a refinement criterion iVj = 10, which insures that the local Jeans length 
is resolved by at least 10 cells. The initial grid contains 256 3 cells and we use 10 levels of 
refinement for an effective resolution of 262144 3 and a minimum grid size of ~ 2.16 AU. 
We apply periodic boundary conditions for hydrodynamics and gravity, and we impose a 
temperature of 10 K for the radiation at the edge of the box. 
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Table 1: Models parameters. 



Model 


n 


«turb 


Ax min (AU) 


Coarse grid 


to (Myr) 


SPHYDRO 


oo 


~ 1(T 5 


2.16 


128 3 


0.4786 


MU130 


~ 136 


~ 0.2 


2.16 


256 3 


0.4935 


MU5 


~ 5.3 


~ 0.2 


2.16 


256 3 


0.5397 


MU2 


~ 2.3 


~ 0.2 


2.16 


256 3 


0.5982 



3. Results 

3.1. Qualitative description 

In this section, we qualitatively describe the results of the four calculations. 
Calculations are synchronized at the time to when the maximum level of refinement is 
reached (see tabled]). The subsequent evolution after t strongly depends on th e physic al 



cond itions (rotation, outflows, etc.). At to, the first hydrostatic core (FHSC, 



Larson 



19691 ) is forming. As expected, we note that the stronger the magnetic field is, the later the 



collapse occurs because magnetic fields "dilutes" gravity. In Figured], column density and 
local Jeans length maps of the four calculations are shown in the xz-plane. The boxes are 
centered at the maximum density of the total computational domain. In the SPHYDRO 
model, only one central fragment is formed and the collapse is nearly spherical. The mass 
of the fragment, i.e. where p > 10 11 cm -3 , at time ~ t + 2.6 kyr is ~ 0.2 M . The 
integrated mass of the envelope from the fragment which is stable against fragmentation, 
i.e. Mi n t/Mj > 1, is ~ 30 M© and the accretion rate is ~ 10~ 4 M© yr _1 . 

In the MU130 model, the additional turbulent support clearly favors fragmentation (6 
fragments are formed at that time) over a region of 2000-4000 AU. The mean separation 
between the fragments is about 1000 AU and they are distributed at this early time along a 
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Fig. 1. — Top: Column density maps integrated in the y-direction for the four models: 
SPHYDRO at time ~ t + 2.6 kyr, MU130 at time ~ t + 25 A kyr, MU5 at time ~ t + 7.5 
kyr, and MU2 at time ~ to + 7.3 kyr. Bottom: Local Jeans length and velocity field cut in 
the xz-pl&ne for the same calculations and at the same time as in the upper row. 

filamentary structure. This separation corresponds to the typical Jeans length associated 
with the region surrounding each fragment. The overall accretion rate on the fragments is 
relatively low, ~ 10~ 5 M Q yr _1 , and consistent with the ones obtained in the low- mass star 
formation framework. The mean mass of the fragments is about 0.2 M Q (~ 1.2 M Q in total) 
and corresponds to the local Jeans mass. 

In the MU5 model, the core has fragmented into two main objects, one being formed 
by the merger of secondary fragments. The fragmentation zone has the same filamentary 
morphology as in the MU130 model, but the extent is smaller. The angular momentum is 
more efficiently transported than in the MU130 model because of the stronger magnetic 
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field (IHennebelle et al.ll201l[ ). The region where rotation can support the collapsing cloud 
is thus smaller. The local Jeans length raises substantially indicating that the core is less 
prone to fragment. The accretion rate in this model is ~ 10 -4 ' 2 M yr~ x . The mass within 
each fragment is ~ 0.25 M©. 

In the MU2 model, there is only one fragment which drives a low velocity outflow 
(^out ~ 2 km s" 1 ) of extent ~ 1500 AU. The accretion rate is ~ 10" 4 ' 1 M Q yr _1 . The local 
Jeans length is much larger than in the MU130 and MU5 models and shows some typical 
features of magnetized dense core collapse. The region of small Jeans length (Lj < 3000 
AU) corresponds to the pseudo-disk that has formed perpendicular to the magnetic field 
lines at larger scales. The feature that appears in the vertical direction corresponds to the 
outflow that has been launched around the central fragment. As radiation preferentially 
escapes perpendicular to the pseudo-disk where the optical depth is smaller, the outflow 
region heats up and the local Jeans length raises. At time ~ t + 7.3 kyr, the mass of 
the fragment is ~ 0.25 M and the associated Jeans stable mass integrated through the 
envelope is ~ 10 M and extends to a radius of ~ 20000 AU. By contrast in the MU5 model, 
the stable Jeans mass associated to each fragment is ~ 1.2 M and extends to a radius of 
~ 800 AU. We note that we have performed lower resolution calculations (Ai m i n ~ 32 AU 
and Aa; m i n ~ 8 AU) which do not show any fragmentation at latter times, while similar 
MU130 simulations show a lot of fragments (see fig. [2]). At times given in fig. |2l the mass 
contained in the fragments with Aa; m i n ~ 32 AU is 7.1 ~ M for the MU130 model and 
~ 1.2 M for MU2 (respectively ~ 0.6 M and ~ 2.2 M with Ax min ~ 8 AU.) 



3.2. Quantitative analysis of the SPHYDRO model 



In this section, we focus on the thermal behavior observed in the SPHYDRO model, 
i.e. without turbulence and magnetic field. Figure |3] shows radial profiles of density and 
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temperature at time y to + 2.6 kyr . The density pr ofile exhibits the classical R 2 slope in 



the envelope (ILarson 



1969 



Penston 



1969 



Shu 



19771 ). Unlike the low mass star formation 



case, for which T oc R~~ 5 in the optically thin envelope, we find that the temperature 
profile in the preshock region ahead of the FHSC is steeper. The optically thick region 
extends up to a radius of ~ 130 AU (vertical grey line) much larger than the FHSC radius 
Rf c ~ 20 AU. The radiation is thus trapped in an optically thick bubble where the infalling 
gas can be efficiently heated and the Jeans length raises efficiently. Compared to the 
low-mass star case, the accretion rate and thus the post-accretion shock temperature are 
larger (~250 K compared to ~70 K). Since /% oc T 2 for temperature < 100 K, this explains 
the larger optical depth found towards higher masses. The temperature profile is well 
fitted in the optically thick region with T oc R~ 7 ^ 8 (dotted line). Such a profile is obtained 
assuming that the rate of change in kinetic energy, ~ GMM/R equals the radiative flux, 
~ R 2 c/(3k,rp)clrT 4 /R (with /% constant for temperature ranging from 100 K to 1000 K). 
At larger radii the optically thin behavior, T oc R~ 05 (dashed line), is recovered after a 
transition regime, until the equilibrium temperature of 10 K is reached. 



3.3. Comparison with the other models 



Figure H] shows the temperature distribution as a function of the density at two times for 
the four models. To visualize the difference between our w ork and other studies, we overplot 



the barotropic law that is used in 



Hennebelle et al. 



(1201 ll ). With a full RMHD model, we 



obtain a spread in temperature and Jeans mass which depends on the magnetization of 
the core that a barotropic EOS cannot reproduce. The SPHYDRO and MU2 models are 
strikingly similar in the heating efficiency, even though the spread in temperature is larger 
in the MU2 model because of the initial turbulence. The rotation in not large enough to 
support the collapsing core at larger radii because of the magnetic braking, and matter 
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Fig. 2. — Low resolution runs for the MU130 and MU2 models. 



is accreted onto a single fragment of similar size and accretion rate as in the SPHYDRO 
model. The accretion luminosity on the FHSC, L acc = GMf c M / R{ c , is thus about the same. 
In the MU130 model, the accretion rate is much lower and the contraction of the FHSC 
slower because of lower temperature. As a consequence, the typical Jeans mass is about 2 
orders of magnitude lower than in the SP HYDRO and MU2 models. As in other models 



( jKrumholz et al. 



2009 



Peters et al. 



2011 



Hennebelle et al. 



20H| ), the MU130 model shows 



fragmentation in several objects even at time to + 6 kyr, where three objects formed with 
separation of 3000-5000 AU. As the magnetic braking is very small here, the radius at which 
the collapse is stopped and where fragmentation is taking place, is thus much larger than 
the FHSC radius, and the accretion luminosity negligible. Even with a fragment undergoing 
second collapse, the radiative feedback would not be efficient at a distance of d ~ 3000 AU 
to suppress fragmentation. For an optically thin envelope, the temperature indeed scales as 
T ~ {Lj (Anad 2 )) 1 ' 4 , where L is the accretion luminosity on the protostar, L = GM+M / R* 
(~ 13K at 3000 AU with M = 1 x 10" 5 M /yr, M* = 0.2 M , and R* ~ 5 R ). Finally, the 
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Fig. 3. — Radial profiles of density and temperature for the SPHYDRO model at the same 
time as in figure [TJ The grey vertical lines indicates the radius at which r = 1 . 

MU5 model shows an intermediate behavior. Although the heating owing to the accretion 
shock on the central fragment is already large, external regions have time to cool and to 
continue collapsing to form a secondary fragment. 
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Fig. 4. — Temperature-density distribution for the four models. The top row corresponds to 
snapshots just before grain evaporation for the SPHYDRO, MU5, and MU2 models, and to 
time to + 6 kyr for the MU130 model. The bottom row corresponds to the same time as in 
figure [TJ The color coding indicates the mass in M Q per bin of equal density and temperature. 
The black lines represent iso- Jeans mass curve, ranging from 10~ 4 M (bottom right line) 
to 10 M (top left dashed line). The yellow curve represents a classical barotropic EOS. 

4. Towards massive star formation? 
4.1. Discussion 

As first suggested by 



Commercon et al. 



(120101 ) in the low mass star formation 



framework, we show that there is a strong interplay between magnetic field and radiative 
transfer which may suppress initial fragmentation. Because of the rotation slowdown owing 
to the magnetic braking (most efficient at scales < 200 — 300 AU), the infall velocity on 
the fragment is then greater. For massive dense cores, the accretion is larger than for low 
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mass ones and the amount of energy radiated away at the FHSC surface, the accretion 
luminosity, is thus much larger. Depending on the initial conditions, we find three totally 
different behaviors: 

• For a spherical massive dense core, without turbulence and magnetic field, we find 
that the core's envelope heats efficiently (SPHYDRO model), because the radius at 
which the accretion luminosity is released is small (FHSC radius). We also find that 
contrary to the low-mass regime, the optically thick regions is much more extended, 
up to ~ 150 AU (~10 AU for the low-mass case). Only one fragment is formed; 

• In the MU130 model, the magnetic braking is too low to transport outwards angular 
momentum produced by the initial turbulence. 7 fragments have formed by the end of 
the calculations. The fragmentation zone extends over a few thousands AU, at which 
accretion luminosity is first released through a radiative shock. Disks with radius 

~ 100 — 200 AU are also formed around the fragments, which give a second accretion 
shock at the disk edges. The infalling gas thus encounters several accretion shocks 
before being accreted onto several FHSCs, which lowers the accretion luminosity. 

• For highly magnetized and turbulent massive dense cores (MU2 model) we find 
very similar results to the SPHYDRO model, i.e. a single fragment, because of the 
magnetic braking that lowers the radius at which the accretion luminosity is released 
to the FHSC radius. 

• We find an intermediate behavior for lower magnetization degrees (MU5 model), 
where the core fragments into two objects of similar properties. The latter model may 
be the most realistic one, in accordance with observations (see hereafter). 
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4.2. Astrophysical consequences 



Our results suggest that the combined effect of magnetic fields and radiative transfer 
could control the early fragmentation of the core and we speculate that this could lead 
either to the formation of isolated massive stars or OB associations. The MU2 model 
gives raise to a single fragment with a large reservoir of mass ~ 10 M Q stable against 
fragmentation. The strong magnetic field case is thus a preferred scenario for non-runaway 
massive stars (i.e. not ejected form a cluster) that are found in isolation. Our proposed 
scenario does not exclude the formation of close massive binaries and OB associations. 
Fragmentation in massive binary system can possibly occur in more massive cores or in cores 
with different initial conditions than ours (for example with stronger density fluctuations 



initially). They could also form during t he second collapse p 



the low mass star formation framework ( Machida et al. 



rase as it has been shown in 



20081 ). OB associations can also 



form by global collapse of a giant molecular cloud containing several massive magnetized 
dense cores. In addition, the MU5 model produce two fragments, that are also associated 



with a relatively 



(IKrumholz et al. 



high 



2007 



eans mass reservoir. Co ntrary to previous lower resolution studies 



2009 



Peters et al. 



201 lh . the secondary fragment is not produced by 



disc fragmentation, but rather by collapse along a filamentary structure. One can expect 



that this early fra gmented system will a. 



following results of 



Bonnell & Bate! (120051 ) . 



so give rise to a close massive binary system 



Recently 



Girart et al. 



(120091 ) report observations of a hot molecular core (HMC) in the 
massive star-forming region G3141 and conclude that the gravitational collapse of the HMC 
is controlled by magnetic field. They also observe a spin-down in the HMC which suggests 
that magnetic braking is acting and removing angular momentum. Last but not least, they 
infer a mass-to-flux over critical mass-to-flux ratio of 2.7 which corroborates our results. At 



later evolution stages, 



Bestenlehner et al. 



(120111 ) observe a very massive star of ~ 150 M Q 
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in apparent isolation from the massive cluster R136, which could be for med from a highly 



magne tized dense core rather than being ejected form the dense cluster. 



Bontemps et al. 



( 120101 ) report imaging of massive dense cores with high angular resolution. They find that 
fragmentation in massive core tends to lead to fewer high-mass fragments inconsistent with 
a pure gravoturbulent fragmentation (which would correspond to the MU130 model). Their 
results support our findings on the enhanced effect of the magnetic braking and radiative 
transfer when mass increases. 



5. Conclusion and prospects 



In this letter, we propose a new mechanism to suppress initial fragmentation of highly 
magnetized massive dense cores and speculate that it can lead to massive star formation. 
Our scenario differs form previo us work, since it do es not invoke stellar mergers as in the 



competitive accretion scenario (IBonnell et al. 



feedback and hig 


1 CO 


(Krumholz et al. 


2009 



Peters et al. 



200 l[) , nor low m ass protostars' radiative 



umn density ( Krumholz fc McKee 



2011 



20081 ). nor disk fragmentation 



. We investigate the early stages of the collapse 
and fragmentation of turbulent, massive, and magnetized dense cores with RMHD 
calculations. We show that the combined effect of magnetic braking and radiative transfer 



suppresses fragmentation in the case of a strong magnetic field a s it has been s 



hown m 



the low-mass star formation framework (jCommercon et al 



2010 



Tomida et al 



2010|). 



Magnetic braking transports angu 
FHSC is thus larger. As shown in 



ar momentum outwards, a nd the accretion rate on the 



Commercon et al. 



( 1201 lbl ). all the infall kinetic energy 



is radiated away at the first core accretion shock and allows greater heating. The interplay 
between magnetic fields and radiative transfer is independent from the initial conditions 
other than magnetic fields strength. In addition, this interplay can only be caught in 
self-consistent RMHD models with thermal (re-)emission and radiative shocks. This effect 
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is expected to become stronger with appropriate physical models 
and second core formation and for protostellar evolution (see e.g. 



or the second collaps e 



Krumholz et al. 



20071 ). 



To reach the final stellar mass, the s ubsequent evolution of the forming prot ostars can 



then follow a disk accretion scenario (IKuiper et al. 



2011 



Seifried et al. 



2011). Further 



Federrath et al 



20 111 ) 



work will imply detailed fragmentation and resolution studies (e.g 
and the introduction of sink particles to follow a longer dynamical range. We also neglect 
i n this work ionization and protostellar feed backs which influence further time evolution 



( IPeters et al. 



2011 



Cunningham et al. 



201ll). 



As a conclusion, we speculate that highly magnetized dense cores are the seed of 
massive stars and good candidates for massive star forming regions. 



Calculations have been performed on the THEO cluster at MPIA and on the JADE 
cluster at CINES. BC thanks Romain Teyssier and Henrik Beuther for useful discussions. 
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